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(yQ [ Abstract 

^ ! 

. The existing experimental data for the deuteron charge radius are discussed. The data 

I of elastic electron scattering are inconsistent with the value obtained in a recent atomic 

physics experiment. Theoretical predictions based on a nonrelativistic description of 
. the deuteron with realistic nucleon-nucleon potentials and with a rather complete set of 

meson-exchange contributions to the charge operator are presented. Corrections arising 
from the quark-gluon substructure of the nucleon are explored in a nonrelativistic quark 
model; the quark-gluon corrections, not accounted for by meson exchange, are small. 



O . 

^ , Our prediction for the deuteron charge radius favors the value of a recent atomic physics 

^ ' experiment. 

PACS number(s): 21.45.+V, 21.10.Ft, 25.30.Bf, 27.10+h, 12.39.Pn 

1. Introduction 

Deuteron properties are quantities of fundamental importance. The hadronic properties are 
used to calibrate parametrizations of the nucleon-nucleon (NN) interaction. The electromag- 
netic (e.m.) properties are crucial tests of the theory of e.m. currents, but simultaneously 
also of the chosen form of the NN interaction. Although there has been a continuous interest 
in all deuteron observables, the observable root-mean-square (rms) charge radius r^h and 
the related but unobservable rms matter radius have recently received special attention 
^, 0, ^ H, 0, @]. This is mainly due to unresolved discrepancies between experimental 
data and theoretical predictions. These discrepancies have even led to speculations 0) §1 
that quark-gluon effects may surface in the deuteron charge radius. We are sceptical that 
these effects can resolve the observed discrepancies. Nevertheless, possible quark-gluon ef- 
fects in the deuteron charge radius deserve a careful investigation. This is the reason why 
this paper studies general short-range corrections of hadronic and of quark-gluon nature for 
the deuteron charge radius, not discussed before. With respect to the history of the subject, 
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Ref. IP gives a comprehensive review of both experimental and theoretical aspects of the 
deuteron radius. 

Sect. 2 summarizes the existing experimental data of elastic electron scattering and of 
an atomic physics experiment relevant for the deuteron charge radius. It does not contain 
novel material. Despite the review of Ref. M the section is important for this paper, since 



it sets the stage for our calculations; it extends our discussion in Ref. ||TOl- All calculations 
of deuteron charge properties require meson-exchange (MEC) corrections for the charge 
operator. Sect. 3 presents our theoretical predictions. Subsect. 3.1 discusses them in a 
hadronic description of the deuteron. Subsect. 3.2 extends this description by taking the 
quark-gluon substructure of nucleons into account. We give our conclusions on the unresolved 
discrepancy between experimental data and theoretical predictions in Sect. 4. 

2. Experimental Data and their Analyses 

The mean square (ms) deuteron charge radius r^;^ is defined as the slope of the deuteron 
charge monopole form factor Fc/i(Q^) oi elastic electron scattering at zero three- momentum 
transfer Q, i.e., 



ch = -Qj7^Fck{Q')\Q^=0. (1) 



d 

dQ 

Equivalently, it can be derived from the slope of the deuteron longitudinal structure func- 
tion A(Q^) at zero momentum transfer, corrected for a small contribution arising from the 
deuteron magnetic moment Hd measured in units of nuclear magnetons fiN = e/2Mp, where 
Mp is the proton mass. 

The latter correction /i^/2Mp = 0.0163 fm^ arises, since the longitudinal structure function 
y4(Q^) is built up from the charge monopole, charge quadrupole and magnetic dipole form 
factors F,h{Q^), Fg{Q^) and F^{Q^) according to A{Q^) = Fi(Q2) + Q7l8 F^{Q^) + 
Q^/6Mp F^(Q^) [|l^; the form factors are considered in the Breit frame in which = 
— g^, q being the four-momentum transfer. The charge quadrupole and the magnetic dipole 
form factors are normalized to the deuteron quadrupole moment in units of fm^ and to the 
deuteron magnetic moment in units of fiN, respectively. 

Characteristic experimental data for A(Q^) at small momentum transfers are displayed in 
Fig. 1(a). They are obtained in elastic electron scattering; the bearing of the experimental 
charge radius, measured by the isotope shift of the Is — 2s transition of atomic hydrogen, 
on A(Q^) is displayed in Fig. 1(b). In Fig. 1, we use the standard form of a high-resolution 
plot which emphasizes experimental errors and possible inconsistencies between experiments. 
We see that v4(Q^) is at low momentum transfers not determined with the desired accuracy. 



Furthermore, the data of Refs. |12| and |13] show opposing trends as already emphasized 



in Ref. |13]; however, the low-momentum transfer data of Ref. |12] have always been 



considered reliable beyond any doubt, whereas those of Ref. |jT3[ have set a new standard 
of experimental accuracy at intermediate momentum transfers. Thus, the discrepancies 
between both data sets create general worries about the low-momentum transfer data, most 
important for extracting the deuteron charge radius Vch- 

Elastic electron scattering is often simultaneously performed on deuteron and proton 
targets. The advantage of such a simultaneous measurement is that the corresponding ratio 
of deuteron and proton charge form factors, FchiQ"^) / G EpiQ,"^) , can be extracted with a 



2 



Experiment 


Teh [fm] Td [fm] 


Tp [fm] rl [fm^] Vch [fm] 


Berard et al. [15 
Akimov et al. [16 
Simon et al. 12 " 
Simon et al. [12|^ 


2.1256(64) 1.9635(45) 
2.098 (26) 1.935 (28) 
2.1159(65) 1.9540(47) 
2.1193(80) 1.9576(68) 


0.805(11) -0.1134(24) 2.0952(60) 
0.817( 8) -0.1170(18) 2.080 (27) 
0.862(12) -0.1170(18) 2.1160(60) 
0.862(12) -0.1170(18) 2.1190(80) 


Schmidt-Kaler et al. [14[ 


2.1303(66) 1.9685(49) 


0.862(12) -0.1130(30) 2.1303(66) 



Table 1: 



Experimental structure radius and charge radius rch of the deuteron. 

and for an 
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Results are given for three elastic electron scattering experiments [jT2|, |T5 
atomic physics experiment []14| . The experimental values [[12|, |15], |16[ for the structure radius 
Td as defined in Eq. (|^) are listed in column 2; the resulting deuteron charge radii using the 
proton and neutron charge radii employed in these references are given in columns 3 to 5. 
In column 1 we derive updated values for the deuteron charge radius rch starting from the 
experimental structure radius r^ of column 2. For this purpose we use the presently accepted 
value for the proton charge radius, i.e., rp 



0.862(12) fm [|T^ but use the same neutron 
charge radius and the standard Darwin-Foldy term, rj-,p = 0.0331 fm^, as in the respective 



references [12, 15, 16 [. Two analyses of the data of Ref. [12| using 3rd and ^''^ 4th order 
polynomials to fit the ratio FchiQ"^) / G EpiQ"^) are quoted. The atomic physics experiment 
gives the radius difference (r^^ — rp)exp = 3.795(19) fm^; the corresponding values for 



the charge and structure radii are obtained from Eq. (3) by using = 0.862(12) fm [ff? 

^2 



-0.113(3) fm2 



18] and rjjp 



0.0331 fm^ 



smaller systematic error. Thus, the difference r'^f^ — r'^, r^ being the ms proton charge radius, 
is experimentally determined with higher accuracy than the deuteron charge radius itself. 
This is the reason why experimentalists prefer to analyze their data for the rms deuteron 
structure radius r^, defined by 

2 _ / 2 2\ 2 2 /o\ 

'^d — Vch ''^pJ ^DF- I'JJ 

The ms structure radius contains as dominant part the experimentally determined difference 
r^^ — Tp. From this one removes (i) the contribution of the neutron finite e.m. size, 
being the neutron ms charge radius, and (ii) the relativistic Darwin-Foldy correction to the 
nonrelativistic one-nucleon charge operator, rj^p = G'|(0)3/(4M^) = 3/(4M^) = 0.0331 
fm^, where (^^(Q^) is the isoscalar charge Sachs form factor of the nucleon and Mn the 
isospin-averaged nucleon mass. In Eq. (||), r^^^, and are observables and rj^p is a 
well-defined number. The rms deuteron structure radius can therefore be considered an 
honest observable quantity of the same standing as the deuteron charge radius rch itself; in 
many experiments it is determined with even higher accuracy. 

The existing experimental values for and rch are listed in Table 1. The quoted values 



refer to the elastic electron scattering experiments of Refs. u2, O, 16[; the experiment of 



Ref. [|T3| whose data are also shown in Fig. 1 was meant to provide information on A(Q ) at 
intermediate momentum transfers; the authors did not extract the charge radius. The most 
recent experimental result on the deuteron charge radius has been obtained in an atomic 



isotope shift measurement and is also given in Table 1. We base the discussion of the 
paper on the results for rch in column 1 of Table 1. We call these results experimental, 
although they cannot be found in the quoted references. The values quoted there are listed 
in columns 2 to 5. We consider the structure radius in column 2 experimentally well- 
determined. However, in Refs. [|T3 and [|T^] the corresponding deuteron charge radii rch, 
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quoted in column 5, were obtained from Eq. (^) using the proton charge radius available 
at that time. We therefore update the charge radius Vch by using in Eq. (3) the currently 
accepted proton charge radius [0, i.e., Vp = 0.862(12) fm, and the values r„ and tdp as 
employed in Refs. |]T2|, |T3], |TB|. We note an unsatisfactory spread of results in column 1 



arising from the analyses of the considered experiments. 

In addition to the charge and structure radii, theoretical papers often discuss the deuteron 
matter radius , 
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^""^ aJo ^^^^ 



oo 



(4) 



It is defined in terms of the deuteron S- and D-state wave functions, i.e., u{r) and w{r) 
respectively. In Eq. r stands for the relative distance between the two nucleons; a factor 

from the volume element is already included in the wave functions. The matter radius 
is theoretically important because it shows the differences between deuteron wave functions, 
obtained from different potential models for the NN interaction, in a direct and pure, though 
integral form. The matter radius is conceptually interesting in the following respects: 

• Due to the factor in Eq. (^) the matter radius weighs heavily the long-range part 
of the deuteron wave function which is dominated by one-pion exchange. All realistic 
parametrizations of the NN interaction have one-pion exchange tails. The correspond- 
ing matter radii typically show a spread of their values which is smaller than 1% as 
demonstrated later on. 

• All realistic parametrizations of the NN interaction account for the deuteron binding 
energy with high precision. Their account of the proton-neutron spin-triplet scattering 
length at and of the asymptotic S'-wave normalization As is far less precise. However, 
one observes |^, |^ that the variations of at and As are linearily related to those 
of the matter radius r^. At the experimental value of the scattering length at = 
5.419(7) fm according to Ref. [l|, [1^ the matter radius becomes = 1.969(4) fm. In 
contrast, the experimental values for the asymptotic normalization As [jl9[ are not as 
precise as those for the triplet scattering length at and are even in conflict with each 
other. Nevertheless, they also suggest a matter radius consistent with = 1.969 fm, 
though with much larger error bars [|I| . Conversely, assuming the empirical correlations 
between at, As, and the matter radius are physically reliable, the experimental value 
for the scattering length at suggests -through the intermediary of the matter radius- 
an asymptotic S'-wave normalization As = 0.885(2) fm~^/^. 

Though we prefer to discuss the model-independent deuteron charge radius Vch, we admit in 
all fairness that the deuteron matter radius Vm is conceptually important and that it may, 
with the help of the well-established empirical linear relations, even serve as a measure of 
consistency in experimental data and in theoretical predictions for characteristic deuteron 
observables. Those observations are the reason why others would like to extract the deuteron 
matter radius also from the experimental electron scattering data. This is clearly impossible 
in a model-independent way, although it has been attempted often. The results of these 
attempts are given in Table 2. 

The difference between the model-dependent ms matter radius and the model-independent 
ms charge and structure radii, r^^ and r^, i.e., 

2 2i2i2i2 i2 i2 / r\ 

^ch = rm + rp + r^ + r^F + ^so + ^[2] . (5) 
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Analysis 


rch [fm] r-m [fm] 




Allen et al. |20| 

McTavish |21| 
Klarsfeld et al. || 
Mustafa 22 1 
Wong P 


2.1095(218) 1.948 (23) 
2.1204( 67) 1.956 ( 5) 
2.1146( 56) 1.950 ( 3) 
2.1189( 52) 1.9547(19) 
2.1150( 52) 1.9502(20) 


0.0000 
0.0150 
0.0135 
0.0135 
0.0143 


Schmidt-Kaler et al. ||14| 
Pachucki et al. 23 


2.1303( 66) 1.9636(49) 
2.1331( 78) 1.9666(66) 


0.0208 
0.0208 



Table 2: Experimental charge radius r^h and matter radius of the deuteron. 

The first five entries in column 2 are the results of various analyses ||l], 21, 22| of 

elastic electron scattering data |jT2|, |l^ aiming at the matter radius r^- The non-nucleonic 
correction r^j used in these analyses is based on vr-meson exchange only; it is listed in column 
3. Our reconstruction of the charge radius is done according to Eq. (^ and given in column 
1. To calculate Vch from the matter radii r.m |1], |20|, |21], |22[ we use Vp = 0.862(12) fm, the 



neutron charge radius of the respective references, the standard Darwin-Foldy contribution. 



DF 



0.0331 fm^, the spin-orbit correction for the Paris potential 



so 



-0.0015 fm^, and 



the non-nucleonic contribution r^j of column 3. In the last two rows the deuteron charge and 
matter radii are calculated from the experimental charge radius difference 



ch 



— r. 



pjexp 



3.795(19) fm2 Q and {r'^ch-rp)exp = 3.807(26) fm^ [H as determined by the isotope shift of 



the Is — 2s transition in atomic hydrogen. The value of Ref. p3] includes higher order QED 
effects. We calculate the corresponding deuteron matter and charge radii using the presently 

-0.113(3) fm2 



5 ^P 



0.862(12) fm fT^] and r. 



accepted values for the nucleonic radii, i.e. 

ig, r|)p = 0.0331 fm^, r|o = -0.0015 fm^, and the two-nucleon contribution r^j = 0.0208 
fm^ of this paper. The latter value is an average of our results for Atmec for the Paris and 
the four Bonn potentials in Table 3. 
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rd = rm + rso + r[2], (6) 
are summarized by two model-dependent corrections. We call these two corrections r|o and 

The first one rgQ is as rj^p in Eq. (|^) a relativistic correction which follows from the 
spin-orbit part of the one-nucleon charge operator, i.e., 

rlo = -6(2Gf,(0) - GUO))Pd^. (7) 

It is model-dependent, since it depends on the deuteron D-state probability Pd] 2Gf^(0) — 
G|(0) = 0.76, Gf^(Q2) being the isoscalar magnetic Sachs form factor of the nucleon. 
For example, for the D-state probability of the Paris potential, Pd = 5.77%, one obtains 
= -0.0015 fm2. 

The second correction r^j is formally defined by Eqs. (5) and (6). As difference between 
observable and unobservable quantities, is model-dependent. Conceptually, the term r^j 
contains all contributions to the charge operator which are considered physically important 
and which are not accounted for by its single-nucleon parts. These corrections are usually 
of two-nucleon nature; this is the reason for the subscript [2]. Thus, r^j has to contain 
two-nucleon meson-exchange corrections r\.jp^. It may contain relativistic boost corrections 
''^boost included by the standard use of nonrelativistic wave functions. It may contain 
A-isobar contributions r^^ in interaction models with explicit A-isobar degrees of freedom 
which yield AA-components in the deuteron wave function and corresponding charge and 
current operators. It may contain additional short-ranged quark-exchange effects t'qec- 
However, this list is not exhaustive given the ingenuity of theoreticians. We therefore write 
symbolically 

^[2] = "^MEC + "^hoost + ^AA + "^QEC + (8) 

Nevertheless, as long as the model-dependent correction r^2] is small compared with the 
deuteron matter radius r.^, it makes sense to apply the corrections t'^q and r^g] to data 
and extract that unobservable quantity from them. The theoretical analyses of Table 
2 usually combine data from several electron scattering experiments. They make the step 
from the structure radius to an "experimental" matter radius according to Eq. (P). 
However, most authors do not simply use Eq. (P), but reanalyse the original experimental 
data with some theoretical bias. The dependence of these analyses on Pn is only slight. 
The non-nucleonic correction r^j is exclusively assumed to arise from one-pion exchange 
effects; the value r^j = r^^^^ = 0.0135 fm^ due to calculations from Ref. ||2^ is often used. 



Compared with the single-experiment analyses of Table 1 we see that different theoretical 
analyses of elastic electron scattering data, usually based on several experiments, yield more 
consistent and with time even converging values. The deuteron matter radius r.^, the main 
result of those analyses, is listed in column 2 of Table 2; the resulting value of about 1.950 
fm is in conflict, by 1%, with the value of 1.969 fm suggested by the experimental scattering 
length at according to the empirical linear relation. This discrepancy is considered serious. 
The analysis |^ of the Saclay data is not included in Table 2. These data do not 



extend to sufficiently low to have by themselves a major impact on the deuteron radius. 
Nevertheless, Refs. f^, ^ conclude that the data set [|l^ favors a larger matter radius than 
the low-momentum transfer data sets. 

In contrast to the theoretical analyses, we like to reconstruct values for the deuteron charge 
radius, listed in column 1. We obtain consistent values around 2.115 fm for it. However, 
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in the last two rows of Table 2 we also give our reinterpretation of the deuteron radius 
resulting from the atomic physics experiment and from an analysis of higher order 
QED corrections for it. Ref. [^] improves some of the corrections of Ref. [Q. A clear 1% 



discrepancy between the deuteron charge radius extracted from elastic electron scattering 
and from the atomic physics experiment surfaces and has been known for some time |TD|, EH]; 



this is a discrepancy between experimental data. This discrepancy remains unexplained and 
is disturbing. In addition and to the distress of theoreticians, no calculation based on realistic 
deuteron wave functions is able to account for the small deuteron charge radius as obtained 
from the elastic electron scattering experiments. This failure, being also of order 1%, has 
been demonstrated before for the deuteron matter radius where the value resulting from 
the electron scattering data according to Table 2 (r^ = 1.950(3) fm) and the value based 
on the experimental scattering length (r^ = 1.969(4) fm) [|l], using the empirical linear 
relation of theoretical nature, are in serious conflict. The discrepancy between experimental 
data and theoretical predictions lead to the speculations of Refs. [|], ^ |[ that quark-gluon 
effects may show up in that observable. 

The theoretical analyses, quoted in Table 2, require an assumption on the non-nucleonic 
correction r^] of the deuteron radius which is model-dependent. Its meson-exchange cor- 



rection has been calculated in the past in |24| for currents of pion-range only and in rather 



rough approximation. Though we expect corrections arising from heavier mesons to be rather 
unimportant, we conclude that a new determination of meson-exchange contributions to rj^j 
is necessary and timely. Furthermore, because of the recent speculations on exotic effects 
we also study the importance of quark-gluon corrections in r^j using a nonrelativistic quark 
model of the deuteron. Both calculations are done in Sect. 3. 

3. Theoretical Predictions for the Deuteron Charge Radius 

The theoretical predictions are displayed in the form 



rch = Vl + rl + rl + r^^ + r|o + r|i]^ (9) 



r.. = Id + r= + aU (10) 

rch = [rl + rl + rl]^ + Atdf + Arcjo + Arp]. (11) 
Eq. (PUP exploits the smallness of the corrections r|,^, rl^ and rL; the equality holds better 



than 0.01% and this accuracy equals the computational accuracy. Eq. (|TT]) defines the 
corrections Arop, A'^so and Ar[2] in an obvious way; these corrections will be recorded in 
this paper. We prefer the linear form of the corrections, since the effect of contributions to 
the deuteron radius can be traced in a transparent way. However, the linear form also has a 
clear conceptual disadvantage: E.g., r'j-)p is a model-independent number whereas Arj^p = 
''^'niFl'^Vrn + '^p + ^n]^ becomes model-dependent due to the denominator. Anyway, in the 
tables below, r'j^p-, r'^so ^^"^ 

always be recovered from the knowledge of ['"m + + '"n] ^ • 



3.1. The Deuteron as a System of Hadrons 

This subsection considers the deuteron as a system of hadrons. Only nucleon degrees of 
freedom are kept active, isobar and meson degrees of freedom are frozen into instantaneous 
potentials and instantaneous one- and two-nucleon currents. Nonrelativistic quantum me- 

,2 

[2] 



chanics is used; a relativistic boost correction is applied. Thus, the nontrivial correction rf, 



is assumed to arise from meson exchange and relativistic boost corrections, i.e. 
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RSC 


Paris 


BonnQ 


BonnA 


BonnB 


BonnC 


TT-pair 


0.0050 


0.0036 


0.0043 


0.0043 


0.0047 


0.0050 


rj-paii 


-0.0001 


-0.0001 


-0.0001 


-0.0001 


-0.0001 


-0.0000 


p-pair 


0.0001 


0.0002 


0.0005 


0.0005 


0.0003 


0.0001 


cu-pair 


-0.0000 


-0.0000 


-0.0001 


-0.0001 


-0.0001 


-0.0001 




0.0004 


0.0003 


0.0001 


0.0001 


0.0002 


0.0002 


TT-ret 


0.0002 


0.0001 


0.0001 


0.0001 


0.0001 


0.0002 


?7-ret 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


p-ret 


0.0000 


0.0000 


0.0001 


0.0000 


0.0000 


0.0000 


cj-ret 


-0.0004 


-0.0004 


-0.0004 


-0.0005 


-0.0004 


-0.0004 


(T-ret 


0.0005 


0.0005 


0.0005 


0.0005 


0.0005 


0.0005 


6-Tet 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0001 


-0.0001 




0.0057 


0.0042 


0.0050 


0.0048 


0.0051 


0.0054 



Table 3: Meson-exchange contributions to the deuteron charge radius. The in- 
dividual contributions are given in fm for six two-nucleon potentials in the normalization 
Ar of Eq. (|l3); they are listed according to the exchange process (pair or ret as defined 
in Ref. [^]) and according to the exchanged meson. The row pTC-f corresponds to the only 
meson-nondiagonal process taken into account. The last row ArMEC contains the sum of all 
contributions. The last digit quoted in the entries is numerically unstable; however, trends 
reflected in the last digit appear correct. 



^"[2] ~ ^MEC + ^ boost- (12) 

The deflnition of meson-exchange currents is based on the extended S-matrix method of 
Ref. and our use of it is described in Ref. p8|. Consistency of the meson-exchange 

currents with the underlying two-nucleon potential in form, i.e., with respect to the hadron 
form factors and with respect to the off-shell extrapolation parameters (/i, z/), and in meson 
content is required. Calculations are carried out for six two-nucleon potentials, i.e., for the 
phenomenological Reid soft-core (RSC) ||2^, for the semiphenomenological Paris [0 and for 
four meson-exchange Bonn [3l]| potentials. The meson-exchange currents used correspond 
to all mesons which the Bonn potentials employ, i.e., the pseudoscalar mesons vr and rj, the 
vector mesons p and uj and the scalar mesons a and 6. The required isoscalar charge operator 
is expanded in powers of (p/M^y), p being a typical nucleon momentum in the nucleus, and 
its contributions up to the relativistic order (p/Mjy)^ are retained. The physics content of 
the meson-diagonal contributions contact/pair (pair) and retardation (ret) is explained in 
Ref. [^]; among the meson-nondiagonal contributions only the pirj part is kept. In this 
paper, we have not calculated the effect of explicit AA wave function components 



on 



the deuteron charge radius. Taking A A components explicitly into account, Ref. sees 
an increase of the deuteron matter radius, Ref. however, a decrease due to subsequent 
meson-exchange corrections; a further study of that cancellation is indicated. 

Our results are collected in Tables 3 and 4. The break-down of the meson-exchange cor- 
rection Atmec into individual contributions is given in Table 3. Their potential dependence 
is moderate, at most 20 %. The dominant contribution arises from the vr-contact term in 
pseudovector vrA^A^ coupling, which is equivalent to the pair term in pseudoscalar coupling. 
Previous calculations often used contributions of pion range only. The result mostly used 
is that of Ref. which is based on the Paris potential. We agree with that result within 
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RSC 


Paris 


BonnQ 


BonnA 


BonnB 


BonnC 




1.9569 


1.9717 


1.9684 


1.9692 


1.9689 


1.9675 


V^m + r2 + r2 


2.1118 


2.1255 


2.1224 


2.1232 


2.1229 


2.1216 




0.0078 


0.0078 


0.0078 


0.0078 


0.0078 


0.0078 


^rso 


-0.0004 


-0.0003 


-0.0003 


-0.0003 


-0.0003 


-0.0003 




0.0057 


0.0042 


0.0050 


0.0048 


0.0051 


0.0054 


Teh 


2.1249 


2.1372 


2.1349 


2.1355 


2.1355 


2.1345 



Table 4: Contributions to the charge radius in the hadronic description of the 
deuteron. The contributions are given in fm for six two-nucleon potentials in the nor- 
malization Ar of Eq. (|l5). The boost correction ^ri^oost is kinematical in the case of the 
deuteron; it is also calculated, but turns out to be zero for all retained digits. Thus, the 
non-nucleonic correction Arpj is identical to Avmec given in Table 3, i.e., Arp] = Avmec- 



10%, but our inclusion of more processes and more mesons makes Avmec larger by about 
15%. 

The results for the deuteron charge radius are given in Table 4. We consider the result 
derived from the RSC potential unreliable mainly because the RSC triplet scattering length 
is unrealistically small. In addition, the construction of a consistent MEC operator for a 



phenomenological potential such as the RSC potential would require some extra effort 
Nevertheless, we list the results of the old-fashioned RSC potential as reference for the 
computational comparison with previous calculations. The results for the other potentials 
show a small spread only; they cluster around the value 2.135 fm. They favor the value 



obtained in the atomic physics experiment as Ref. p6[ also concludes for its calculations. 
According to the results of Refs. |3^ this conclusion will remain firm even if effects of 
AA-components in the deuteron wave function were taken into account. 

3.2. The Deuteron as a Six-Quark System 

This section considers the deuteron as a system of six nonrelativistic quarks. The six-quark 
deuteron wave function is fully antisymmetrized; it is obtained according to the Resonating 
Group Method (RGM); only the channel with two asymptotic nucleons is retained. However, 
even a single channel calculation leads to a deuteron wave function which has AA and hidden 
color admixtures in the short-range region M^ . This is a consequence of the Pauli principle 



at the quark level. 

The calculation is based on a chiral-invariant quark hamiltonian with instantaneous two- 
quark potentials; the interaction contains a quadratic confinement potential and residual 
interactions due to gluon-, vr- and cr-exchange. 

Fig. 2 displays the quark-quark interaction; the antisymmetrizer for the six-quark wave 
function is considered as part of the displayed interaction operator. The antisymmetrizer 
acts on the product wave function of two three-quark clusters with nucleon quantum numbers 
correlated by a relative cluster-cluster wave function which is not yet antisymmetrized with 
respect to quarks belonging to different clusters. 

Two models are investigated for the deuteron. In the first one, called model A, the 
quark hamiltonian is taken into account in full within the RGM calculation of Ref. \ ^% . 
Compared with the traditional purely nucleonic two-nucleon potentials, model A yields an 
attractive nonlocal effective two-nucleon interaction at small relative distances which arises 
from (j-exchange with simultaneous quark interchange shown in Figs. 2(b-c) and 2(e-g). The 
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second one, called model B, simplifies the dynamics by retaining only the direct cx-exchange 
between three-quark nucleon clusters shown in Fig. 2(d). In model B the intracluster 
(T-exchange diagram of Fig. 2(a) and the a-exchange diagrams with simultaneous quark 
interchange of Figs. 2(b-c) and 2(e-g) are omitted; the attraction at small relative distances 
is less pronounced compared with model A. 

The comparison of the results of models A and B shows the importance of the attractive 
nonlocality due to the cr-quark exchange diagrams of Figs. 2(b-c) and 2(e-g) for the deuteron 
charge radius. Due to the attractive nonlocality at small distances, the deuteron wave 
function of model A, i.e., the relative wave function between nucleonic three-quark clusters, 
properly renormalized to account for the norm kernel according to Ref. [^, shows an 



increased probability at small relative distances with respect to model B. Compared with 
the purely nucleonic Paris potential wave function, both wave functions of model A and B 
have a slightly increased probability at small internucleon distances but are otherwise similar 
to that of the Paris potential. The six-quark deuteron wave function does not exhibit any 
exotic structures at small relative distances in contrast to the speculations of Refs. 0, |^. 

The charge and current of the deuteron is carried by the quarks. Sample processes 
for the charge operators are displayed in Fig. 3. In addition to the single-quark charge 
operator of Figs. 3(a-c) there are exchange corrections of two-quark nature; they arise from 
gluon-, TT- and a-exchange according to Figs. 3(d-j). Only gluon-and vr-exchange currents 
are considered in the calculation. Currents arising from a-exchange are not taken into 
account; they are assumed to be smaller than those of pion range according to the general 
experience with cx-meson exchange currents between hadrons. The quark processes of Fig. 
3 can be resummed (i) into those of single-nucleon character, (ii) into those of two-nucleon 
TT-exchange character, and (iii) into novel ones of quark-exchange character, called quark- 
exchange currents (QEC). Especially, the relative weight of these three processes is important 
for the discussion of this paper. The quark-exchange current is defined as the sum of all 
processes calculated with the proper six-quark wave functions minus the processes (a), (d) 
and (g) calculated with renormalized wave functions. In other words, the difference between 
processes (a), (d) and (g) calculated with unrenormalized and renormalized six-quark wave 
functions has to be added to the contributions of diagrams (b), (c), (f), (h-j) in order to 



obtain the proper quark exchange current |^ . A more detailed discussion of the charge and 
current model is given in Ref. . 

Our results are given in Tables 5 and 6. Table 5 gives a break-down of the quark-exchange 
processes into individual contributions. Table 6 sums the results up to the charge radius. 
The difference in the predictions for models A and B shows that the short-ranged nonlocal 
attraction due to Figs. 2(b-c) and 2(e-g) indeed leads to some redistribution of matter 
towards shorter distances which affects already the deuteron matter radius and therefore the 
impulse approximation to the charge radius. We find that the novel quark-exchange current 
contributions to the deuteron charge radius amount to about 20% of the traditional tt- 
meson exchange contribution; they are therefore of the same size as the non-pionic exchange 
current corrections considered in Subsect. 3.1. The e.m. QEC corrections are, however, 
smaller than the quark exchange effects generating the attractive nonlocality responsible for 
the redistribution of matter to smaller distances. 

With respect to the speculations that quark-gluon effects may be important for the the- 
oretical prediction of the deuteron charge radius, we therefore have to conclude that our 
calculations do not support them. As expected, the hadronic description of the deuteron 
together with the meson-exchange corrections of Subsect. 3.1 remains valid. Furthermore, 
the discrepancy between the deuteron charge radius as extracted from elastic electron scat- 
tering data and its theoretical prediction remains also in a quark model description of the 
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Model A Model B 


impulse 


0.0001 0.0001 


TT-pair 


0.0003 0.0004 


gluon 


0.0004 0.0004 




0.0008 0.0009 



Table 5: Quark-exchange contributions to the deuteron charge radius. The con- 
tributions are given in fm for the two quark models A and B in the normalization Ar of Eq. 
(P^ip. Row 1: impulse charge operator with quark exchange according to Figs. 3(b-c). Row 
2: pion-pair charge operator with quark exchange according to Figs. 3(e-f), 3(h-j). Row 
3: gluon-pair charge operator with quark exchange according to Figs. 3(e-f), 3(h-j). The 
last row /S.rQEc contains the sum of all contributions. The last digit quoted in the entries is 
numerically unstable; however, trends reflected in the last digit appear correct. 





Model A Model B 




1.9657 1.9680 


\/rl + r^ + rl 


2.1199 2.1220 


Atdf 


0.0078 0.0078 


^rso 


-0.0004 - 0.0004 


Atmec 


0.0043 0.0044 


Atqec 


0.0008 0.0009 


Teh 


2.1324 2.1348 



Table 6: Contributions to the charge radius in the six-quark description of 
the deuteron. The contributions are given in fm for the two models A and B in the 
normalization Ar of Eq. (p!TD; they are graphically displayed in Fig. 3. The processes of 
Figs. 3(a) and 3(d) correspond to the conventional ones of a single-nucleon charge operator; 
the results are given in the first four rows; the experimental values Vp = 0.862(12) fm |p!7| 
and = —0.113(3) fm are used. With the small quark core radius b = 0.5184 fm used 
here the experimental proton and neutron charge radii are not exactly reproduced by the 
underlying quark model; see however ||3^. The meson-exchange correction ArMEC accounts 
for the process of Fig. 3(g); it corresponds to the canonical pion-pair term in the hadronic 
description of the deuteron. The processes of Figs. 3(b-c), 3(e-f), 3(h-j) have no counterpart 
in the hadronic description; they yield the novel quark exchange contributions Atqec] their 
individual contributions are listed in Table 5. 
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deuteron. 



4. Conclusion 



We recall: For a long time the electron scattering data of Ref. |T^, yielding a deuteron 
charge radius of about 2.115 fm, were considered valid beyond any reasonable doubt. In 
fact, the value for the proton charge radius [|^, still unquestioned as standard, is derived 



from the same experiment. Thus, the inability of the purely nucleonic deuteron description 
to account for a deuteron charge radius of 2.115 fm - the theoretical predictions clustering 
around 2.135 fm are off by 1% - was taken to be a pure theory problem. Speculations on 
exotic effects surfacing in the deuteron charge radius started to flourish, though the deuteron 
size is dominated by the wave function tail, the regime in which the two nucleons are well 



separated. The appearance of the electron scattering data of Ref. [12] changed that situation. 



In principle, the experiment of Ref. [|13| cannot provide very detailed information on the 
deuteron charge radius; it measures the structure function v4(Q^) at intermediate momentum 
transfers, but there it is in disagreement with Ref. [1^; thus, the suspicion arose that the 
discrepancy between experimental value and theoretical prediction may be rooted in an 
experimental problem after all. This suspicion gets now further support by the result of the 



atomic physics experiment |TJ] which suggests a larger experimental deuteron charge radius 
of 2.130 fm consistent with the theoretical value. When making this conclusion we assume 
that e.m. corrections beyond one-photon exchange are - on the level of a 1% comparison - 
the same for a bound and a scattered electron; in fact, this holds true in case of the proton 
for which electron scattering and an atomic physics experiment |^ yield consistent data. 
The recalled historic evolution of the problem provides the background for the present paper. 

This paper presents a recalculation of the deuteron charge radius within a purely nucleonic 
picture of the nucleus, with a complete set of meson-exchange corrections included in the 
charge operator. A corresponding calculation within a six-quark picture of the deuteron 
yields only a tiny additional correction of the order of 0.1% arising from processes of true 
quark nature, which are not accounted for in the hadronic description. The size of this 
quark correction is consistent with a previous rough estimate . As expected, the hadronic 
description of the deuteron radius is valid. Quark-gluon effects are very small for the static 
charge properties of a nuclear system bound as lightly as the deuteron. 

Our theoretical predictions are based on six realistic two-nucleon potentials and on a 
nonrelativistic constituent quark model. All yield a value of about 2.135 fm for the deuteron 
charge radius. Only the value for RSC is, with 2.125 fm, significantly smaller; we mistrust 
the RSC value at this level of precision for the reasons discussed in sect. 3.1. Our theo- 
retical predictions of about 2.135 fm disagree with the experimental values of around 2.115 
fm obtained in the analyses of elastic electron scattering data and quoted in Table 2. We 
therefore also present our predictions for the deuteron longitudinal structure function A(Q^) 
at finite momentum transfers Q in Fig. 4 which shows how the two principal data sets []12| 
and are in conflict. This conflict suggests that a remeasurement of ^(Q^) in elastic 



electron scattering, especially at low momentum transfers, is necessary. According to Fig. 4, 
our predictions can describe the data set of Ref. [|T^] which, however, has only little impact 
on the deuteron charge radius, and therefore was left out from Table 2. Nevertheless, the 



data of Ref. |]13l have a trend that favors a larger charge radius consistent with the atomic 
physics experiment as was already noticed in Ref. |25|. This conclusion is derived from 
the observation that all theoretical predictions for A(Q^) in the intermediate momentum 
transfer range of Fig. 4 are ordered according to their values for the charge radius: Exper- 
imentally, data in the intermediate momentum transfer regime provide no information on 
the charge radius; theoretically, however, the underlying hadronic force and current model 
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seems to provide a rather rigid relation between the e.m. properties at low and intermediate 
momentum transfers. 

Finally, we note that the theoretical prediction of this paper for which includes a 
complete set of exchange current operators and quark-gluon corrections is in agreement with 
the deuteron charge radius obtained in the recent atomic physics measurement 



Added Note 

After submitting this paper, a reanalysis of the relevant data on elastic electron-deuteron 
scattering including two-photon exchange corrections appeared The newly extracted 



deuteron charge radius rch = 2.128(11) fm |^ is consistent with the atomic physics experi- 
ment and theoretical predictions. 
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Figure 1: (a) Experimental data for the longitudinal deuteron structure function y4(Q^). 
A high-resolution representation Ay4(Q^) = (y4(Q^) — AjA,parisiQ'^)) / AiA,ParisiQ,'^) * 100 is 
chosen. AjA ParisiQ,'^) is the theoretical result derived from the Paris potential in nonrela- 
tivistic impulse approximation using the nucleonic Dirac form factor (Q^) instead of the 

(Q^) in the nonrelativistic charge operator. The data are plotted in 
— g^, q being the four- momentum transfer. The two 



Sachs form factor 
the Breit frame function of 



elastic electron scattering data sets of Refs. [|T^ (bullets) and (diamonds) appear to be 
inconsistent; for this conclusion data up to = 5 fm~^ are shown. 

(b) An attempt is made to indicate the importance of the deuteron charge radius for the de- 
termination of the structure function A(Q^) at small momentum transfers. The lower dashed- 
dotted curve refers to the linear approximation v4(Q^) = (1 — r^^^Q^/G)^ with Vch = 2.1303 
fm from the atomic physics experiment . The upper dashed-dotted curve shows the cor- 
responding linear approximation for the charge radius = 2.1146 fm resulting from the 
analysis of the elastic electron scattering experiments and |TB[. The reference result 
AiA,paris{Q'^) is linearized in the same way for these theoretical curves. We note that the 
charge radius 2.1303 fm of Ref. [|14| is inconsistent with the low momentum transfer data 
of Ref. [1^, but seems to be consistent with the Saclay data. In the context of Fig. 4 we 
shall note that the linear approximation is poor at momentum transfers larger than 0.5 
fm~^. Thus, the extrapolation of the elastic electron scattering data to the slope of A(Q^) 
or -Fch(Q^) at zero momentum transfer has to be done with great care. 
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Figure 2: Sample contributions to the two-nucleon interaction arising from meson-exchange 
in a quark model description of the deuteron. The horizontal dashed line stands for the 
meson; cx-exchange is chosen here as an example. Process (a) is the intracluster term, 
process (d) the direct intercluster term and the other processes are quark-exchange terms. 
The quark model A includes all a-exchange processes, quark model B only the intercluster 
cr-exchange process (d) without quark exchange terms. In contrast to Ref. ||3^ both models 
A and B employ cr-exchange between quarks and not between the centers of three-quark 
clusters. 
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Figure 3: Sample contributions to the one-nucleon and two-nucleon currents in a quark 
description of the deuteron. The current has single-quark and irreducible two-quark parts. 
The single-quark current can be split into direct and quark-exchange pieces according to 
processes (a) and (b), (c), respectively. The two-quark current arises from tt- and gluon- 
exchange; only the 7r-exchange is indicated as horizontal dashed line in the figure. The 
two-quark current can be split into direct intracluster (d), direct intercluster (g) and the 
quark-exchange processes (e), (f), (h)-(j). Processes (a) and (d) correspond to the single- 
nucleon current and process (g) to the two-nucleon pion-exchange current in a hadronic 
description when calculated with renormalized six-quark wave functions. The break-down of 
the quark-exchange processes into quark-exchange impulse, quark-exchange pion and quark- 
exchange gluon given in Table 5 is now obvious. 
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Figure 4: Longitudinal deuteron structure function y4(Q^). A high resolution representation 

AA(Q2) = (A(Q2) - AiA,ParisiQ^))/AiA,Paris{Q^) * 100 is choseu. AjA,ParisiQ'^) IS the 

theoretical result derived from the Paris potential in nonrelativistic impulse approximation 
using the nucleonic Dirac form factor F^^(Q^) instead of the Sachs form factor G'|(Q^) for 
the nonrelativistic charge operator. Theoretical predictions for three realistic two-nucleon 
potentials RSC, Paris and Bonn C, shown as dotted, solid, and dashed curves, respectively, 
are compared with the data. The results for the Bonn A and B potentials cluster around the 
Bonn C curve; our prediction for the charge radius of the Bonn C potential, = 2.1345 
fm (see table 4) is consistent with the atomic physics experiment (table 2). We note that 
the linear approximation for Ay4(Q^), used in Fig. 1 to emphasize the different trend for 
AAIQ"^) arising from different charge radii, is not valid anymore for momentum transfers 
larger than 0.5 fm~^. We observe that the theoretical predictions for the Paris and for all 
Bonn potentials cannot explain the data of Ref. however, they appear to be consistent 
with the data set of Ref. . 
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